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Abstract 

We constrain the masses of scalar-tops (stop) by analyzing the new precision Teva- 
tron measurement of the VU-boson mass and the LHC/Tevatron indications of a 
Higgs boson of mass 125. 5±1 GeV. Our study adopts Natural SUSY with low fine- 
tuning, which has multi-TeV first and second generation squarks and a light Hig- 
gsino mixing parameter u =150 GeV. An effective Lagrangian calculation is made of 
m/j to 3-loops using the H3m program with weak scale SUSY parameters obtained 
from RGE evolution from the GUT scale in the Natural SUSY scenario. The SUSY 
radiative corrections to the Higgs mass imply maximal off-diagonal elements of the 
(~^. ■ stop mass-matrix and a mass splitting of the two stops larger than 400 GeV. 
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Supersymmetry(SUSY) is a theoretically attractive extension of the Standard 
Model(SM) that may explain the hierarchy of the weak scale and the Planck 
scale. Of the SUSY particles, the lighter scalar top squark may have a sub-TeV 
mass and be detectable by LHC experiments. Existence of a light top-squark is 
particularly suggested by the Natural SUSY model [lT2T3Tll5ll6ll^^ 
that has less fine tuning. The first and second generation squarks have multi- 
TeV masses to mitigate unwanted flavor changing neutral currents (FCNC) 
and large CP violation. For a third generation scalar GUT-scale mass mo (3) < 
1 TeV, mfj is less than 400 GeV from the running of the RGE equations |17j . 

A light top squark can give a significant radiative contribution to the VU-boson 
mass. The precision of My/ has been improved by recent Tevatron measure- 
ments; M w = 80, 387 ±12(stat.) ±15 (syst.) MeV by the CDF collaboration [22] 
and M w = 80, 3Q7±13(stat.) ±22(syst.) MeV by the DO collaboration [23]. In- 
cluding these measurements, the world average Mw is shifted downward from 
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Table 1 

Uncertainty of the SM Mw prediction from the uncertainties of the parameters. 
Beside these errors, there is another uncertainty due to missing higher order cor- 
rections, which is estimated as about 4 MeV. [27] 

HP M^ p = 80, 399 ± 26 MeV to 80, 385 ± 15 MeV. The SM prediction [25J2S] 
of Mw at 2-loop order is 

M^ M = 80,361 ±7 MeV. (1) 

where we have used the numerical formula of ref.[27] with central values of 
parameters [28] The uncertainties of the SM prediction of Mw resulting from 
the uncertainties of these input parameters are summarized in Table [TJ 

The LHC experiments have reported indications of a Higgs boson at mass 
125.3±0.4 stat ±0.5 syst GeV in CMS dataj2H] and at 126.0 ±0.4 stat ±0.4 syst GeV 
in ATLAS data|30j. Accordingly, we assume a Higgs boson mass of 125.5 ± 1. 
GeV in our study. Then, the difference of the experimental and SM values of 
M w is 

M^ p - = 24 ± 15 MeV. (2) 

As can be seen in Table HJ the largest source uncertainty in Mjy 1 (of 6.0 MeV) 
is from the uncertainty 5m t = 1 GeV in the top mass measurement. It is 
significantly smaller than the experimental uncertainty in M^ p (of 15 MeV), 
given in Eq. (J2J). 

The contributions of SUSY particles to the one- loop calculation of M w [3T] 
along with the W self-energy at the two loop level |32] can account for the 1.6a 
deviation of the experimental value from the SM prediction|31J. Conversely, 
the Mw measurement gives a constraint on the squark masses of the third 
generation, m^, m^ 2 , and mr . We assume no mixing in sbottom sector since 
that off-diagonal element is proportional to m;>; is irrelevant to 8Mw- 

The dominant SUSY radiative corrections to m h are due to loops of t\ and t 2 - 
Implications of a 125 GeV Higgs boson for supersymmetric models are investi- 
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gated in ref . |33j . If rrih is confirmed with the value of the present Higgs-boson 
signal ~125.5 GeV, the values of , mf a and the top squark mixing angle Of 
can be constrained from the measured m^. We investigate how a Higgs mass 
rrih = 125.5 ± 1.0 GeV and the new experimental value of Myy constrain the 
third generation SUSY scalar top masses. 



Constraint from My/ The Myy prediction is obtained by calculating the 
muon lifetime [25. 26, 3 lj. The SUSY correction Ar to the Fermi constant is 



G e 2 

^L= „ „ (1 + Ar) (3) 

where s\y = sm6yy and 6yy is weak mixing angle which is defined by the 
experimental values of W/Z pole mass M w / Z as 



fz 



: 2 w ^cos 2 6 w = 1 ^. (4) 



Ar is calculated [S] in the MSSM, and the corresponding Mw prediction is 
obtained by iterative solution of the equation 
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2 ^4 y/2G„M t 
Then, the correction to M^ is at one-loop level is 



TTCt 

[l+Ar(M w ,M z ,m t , ■■■))}. (5) 



SM^ = -M Z 2 W w 2 Ar . (6) 



c 



w °w 



Ar is given by |25(26] 



c 2 w (SMl 5M 2 \ / a \ 

The first term on the left-hand-side is the on-shell self-energy correction to 
gauge boson masses; -rrf- — rrsF = ^ ■ ^ a 1S the radiative 

M z w z w 

correction to the fine structure constant a . The remainder term (Ar) rem . 
includes vertex corrections and box diagrams at one loop level which give 
subleading contributions compared with the 1st term of Eq. (J7])[31j. 

The main contribution to SMyy is the on-shell gauge-boson self energy, which 
is well approximated [32, 34J with its value at zero momenta as 
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Mfv J = s 2 w 



-^Ap (8) 



where Ap is the deviation of the p parameter due to new physics in the EW 
precision measurements. It is related to the T parameter [55] by 

Ap~a(M z )T . (9) 



The squark, slepton, and neutralino/chargino loops contribute to Ap at 1- 
loop level, which we denote as Apo- The neutralino/chargino contributions 
are small|36j, and the slepton contributions are suppressed relative to squark 
contributions by color, and thus the squark contributions are dominant. It 
is well known [57] that the weak SU(2)i isospin violation from SUSY doublet 
masses gives non-zero contributions to 5Mw- The scalar-top sector is expected 
to have a large L — R mixing since the off-diagonal elements of the top squark 
mass matrix are proportional to m t . Finally, 5Mw is given by |32|l34"] 

M w c 2 w 

Ap0 = i^^ s i c i F ^ m h m2 + c t ?f »K.<) + ^o(m|,mfj](10) 

where F (a, b) = a + b— |z|ln| . = sin^, ci = cos%, and 6f is the top squark 
mixing angle. The 2-loop gluon/gluino exchange effects, Apf^^ n / gluino , are 

neglected since they are subleading compared with the 1-loop Ap for M susy ~ 
300 GeV[34j. The prediction of M w in SUSY is then M w = M^ M + 5M W . 
From Eq. ( ITU]) the 5M W of Eq. (0) corresponds to 



Ap = (4.2±2.7) x 10" 4 , T = 0.054 ± 0.034 (II) 

The uncertainty is substantially reduced from that of the previous global elec- 
troweak precision analyses: Ap = (3.67 ±8.82) x 10" 4 [38], T = 0.03 ±0.11 [39]. 

By using Eq. (TTUj) with (J2J), we can determine the allowed region in the 
mi , Ami plane for a given value of % Here Am ( - = (m^ 2 — raj ). The case 
&i = | is shown in Fig. [TJ Note that X t and % are independent because the 
soft-SUSY parameters in the diagonal elements are different. 

We also note that mg £ in Eq. ( ITU]) is given by , m 4 ~ 2 , and Of 

m 2 = m\cos 2 9- t + m|sin 2 6> t ~ - m\ + m\ - M^cos2/3 . (12) 
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Fig. 1 . Allowed regions in the {mi , Am f -) plane for % = f ; Ami = ( m i 2 ~ m t± ) • Black 
(Red) solid lines are 5M\y = 24 MeV (maximum m^ with X t pea k = — \^6M Bmy ). 
The Blue (dark-shaded) region is m^ = 123.5 to 127.5 GeV and the white line 
represents its central value m^ = 125.5 GeV. The Green (medium-shaded) region 
is allowed by 5M\y at 90% CL, and the dotdashed lines represent its la deviation, 
5M W = 24 ± 15 MeV. 

Equation (12) is symmetric under the exchange 

m h m h , Ct s b i- e - h ->• tt/2 - 9- t . (13) 

Constraint from m^o The mass of the Higgs boson in the MSSM receives 
substantial radiative corrections to the tree level result. The scalar-top sec- 
tor gives the dominant contribution, for which Am\ oc mj/v 2 . Tremendous 
efforts [30]- [71] have been expended to calculate m h with sufficient accuracy 
to compare with LHC measurements, and the Higgs mass has been calcu- 
lated through the 3-loop level, a t a 2 , for the leading (m t ) 4 corrections [MfPo] 
and partially at 4-loop level[59]. The dominant contributions arise from su- 
persymmetric loops involving the top squarks, along with gluon and gluino 
exchanges. 

There are several different approaches that have been used in the theoretical 
evaluation of nih : perturbative calculation of the Higgs self energy diagrams 
to (i) 2-loop and (ii) 3-loop orders, (iii) effective field theory (EFT) methods 
based on second derivatives of an effective Higgs potential, (iv) effective poten- 
tial method based on RGE evolution from the GUT scale, and (v) the effective 
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Lagrangian method. We succinctly summarize the five methodologies: 

i) The FeynHiggs package [58] calculates rrih diagrammatically in 2-loop order 
in the on-shell(OS) renormalization scheme. 

ii) A MATHEMATICA program, H3m[65j, does the three-loop calculation; it 
is interfaced with the 2-loop FeynHiggs program for rrih predictions. A numer- 
ical 3-loop accuracy on rrih has been estimated to be < 1 GeV. However, its 
expansion in mass-squared ratios does not apply in some parameter regions 
relevant to Natural SUSY. 

iii) In the EFT 2-loop leading-log approximation [48 ,43 .67] . m\ is calculated 
in the limit of stop matrix elements M L = M B [671l6^1l6W0lr7T] . and M L > 
M R ^Q\. 

The ml formula in the general case with Mi ^ Mr, is given in the large ttia 
limit by[67j 



m h,EFT2( m t 1 , rn h , x t ) = M 2 z c 2 w + 



2n z v' L 12 

I /3m 2 \ / t 2 + t 2 

+^2\—r- 327ra s (m t ) J x ( X t t max + max min + {2t - t max - t min )t E 



mf mf 



where v = l/yy/2Gp — 246 GeV and the contribution from the sbottom 
sector can be omitted so long as tan (3 is not close to its upper bound of ~ 60. 

In the above equation, X t is related with the stop-mixing parameter X t = 
A t — /xcot/3 by 



X t = 2\Xt\ 



,ln(m| 2 /m| 



2 2 

m 4 - — m~ 



'2 



ln(m 2 lm 2 r % 



*2 



2 



(15) 



In Eqs. (fT4l) and ffT5l) the X t is a quantity regularized with the renormalization 
scale fi = M susy in the MS scheme, while the running top quark mass fh t 
is evaluated at [A — m 4 itself in the MS scheme. m t (fj,) was calculated in 
DR scheme by ref.[72] and in Q(a£) [731174] . Its value in the MS scheme is 
rh t = 163.71 ±0.95 GeV|39] which corresponds to the on-shell top quark mass 
M t = 173.4 ±1.0 GeV. 

The X t in Eq. (IT51) is well approximated as 
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X t = 2x 2 t -^, x t = ^- (10) 
' 6' M susy 



with the choice of SUSY breaking scale 



susy 



(17) 



The m1 EFT2 of Eq. (fl4"|) has its maximum at \x t \ = \(x t ) max 

| = v 7 ^ or \X t \ = 

\(Xt)max\ — y/QM susy , for which X t = 6. It is also a common feature of the 
analytic EFT formula at 1- and 2-loop levels [67.68 69 . 70j . A region \X t \ ~ 
V^Afsusy is theoretically not allowed from considerations of false vacuum of 
charge and color symmetry breaking |75|76f77|Tf8] . 

-^max,min are related to the stop squared-mass matrix M~ in on-shell(OS) 
renormalization scheme as 



M t X° s Ml J \-(ml-ml)c- t s- t m\s\ + m\c\ 



(M OS )" in = max, min{Mi M|} = ta - * ± 



\ 




(M,X t os )2. (19) 



Our sign convention of X t agrees with that used in ref . [70] . X; OS is the on- 
shell stop mass matrix parameter. The relation between M^ y and X® s in OS 
scheme and those in MS scheme are given in [70] , see also [57] . Here we treat 
M°f xmin as being equal to M maXjmin in Eq. (JHJ) since the difference is small 
(less than 4%) for M susy > 1 TeV. 

In Eq. ( ITT)) , the r.h.s is given by the on-shell stop masses and thus, more 
precisely Eq. (117)1 is . Here we regard M^ y as being equal to M susy in 
MS scheme since the difference is small. 



On the other hand, X t affects a relatively large difference between DR and 
OS schemes. Numerically, we define the ratio 



(X t ) meLX /(X° s ) 



(20) 



which is about 1.2 from the formula relating MS and OS schemes give 111 
Carena et al.JTD]. Coincidentally, k ~ a/6/2.0. We choose this form because 



MS 



susy 



-"-'susy 



+ 



+ 



3X t 



where the renormalization prescription is not specified in 0(a & 



in 1-loop level [70 
i term. 
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the factor a/6 matches the Xt value in the MS scheme giving maximum X t of 
Eq. ( Fl6|) which leads to maximum m\ eft2 °f Eq. (I14p . The 2.0 in the denomi- 
nator is given as a numerical value of the ratio (X° s ) max /M susy in ref.JTD]. We 
have also checked the ratio (I2"0~j) by using Isajet 7.83JHD]: Isajet adopts the DR 
scheme and DR ~ MS and converts to OS stop masses using[57]. Isajet out- 
puts of Xf R and on-shell stop masses numerically consistent with the relation 
(Xp*) max /pff s ) max = v / 6/2.0. (See, also, the caption of Figffl) We apply this 
relation (|20|) in the region close to "maximal mixing", |X t |/M susy ~ a/6: 

X t = nX? s , k = a/6/2.0 (21) 

The EFT method is not gauge-fixing invariant [59J. Nonetheless, it is found to 
give a good approximation when compared to other methods. The formulae 
tHM with (fl5l) gives larger m^ values by about 1 GeV than the results of H3m 
with the inputs of the natural SUSY benchmark points, as will be commented 
on below. 

The m\ formula obtained from the 2-loop diagrammatic approach (i) can be 
matched to the EFT formula above by adjusting the renormalization prescription [70J, 
except for additional non-logarithmic terms in the diagrammatic formula that 
give asymmetric heights of the peak m^ at X t > and X t < 0. The latter 
contributions arise from SUSY threshold effects that are not taken into ac- 
count in the RGE running down from the SUSY-breaking scale that includes 
logarithms of M susy /m t . 

iv) In the unification approach, RGEs are evolved from the GUT coupling 
unification scale[72j, where the 1st and 2nd generation scalars in Natural 
SUSY have a mo ~ 10 TeV mass and the 3rd-generation scalars have mo ~ 1 
TeV [TT]r79"] . The Higgs potential at the SUSY breaking scale M susy is based on 
one-loop MSSM radiative corrections that are RGE improved. With the choice 
of M susy = y/m^ mi 2 , the most important two-loop effects [66J are included in 
the effective potential. The RGE evolution is implemented with the ISASUSY 
package [80, 8 1J, with a scan over GUT scale parameters. 

v) In the effective Lagrangian approach, the gauge couplings, the Yukawa cou- 
plings, and the soft-SUSY terms are also RGE evolved to the weak scale from 
high scale boundary values, where the gauge couplings unify. The ISASUSY 
program for this RGE evolution incorporates SUSY threshold effects. [80|8"T] . 
The weak scale parameters so obtained are taken as input to the diagram- 
matic calculation at 2-loop order by the FeynHiggs[58j or 3-loop order by 
the H3m[65j. It has been argued [59] that this method may provide the most 
accurate evaluation of the leading and next-to-leading contributions to m^ 
in 3-loop order in the approximation of large QCD and top-quark Yukawa 
couplings. 
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We adopt the latter approach in the framework of natural SUSY using IS ASUS Y [80] 
with a scan over GUT scale input parameters. We have also corrected the sign 
convention of X t in ISASUSY in order to match ours. We then evaluate 
using the H3m program with the ISASUSY input for the SUSY parameters at 
the weak scale. Specifically, we adopt the benchmark line NS3 of Ref. [17] that 
has a Higgsino mass term \i = 150 GeV and other Natural SUSY benchmark 
points RNS1 and RNS2 of Ref pj. 0- The NS3 gives m h = 123.5 GeV that 
is consistent with the LHC experimental value. There is a strong preference 
for A t (M susy ) > and tan/3 > 10 in Natural SUSY[17]. Since /i is small in 




Fig. 2. At(M susy ) dependence of rrih in 3-loop calculation by H3m with the effec- 
tive Lagrangian method. (Solid circles). The input parameters are a natural SUSY 
benchmark line (NS3): (mj B> m t2 b) = (812.5,1623.2) GeV which corresponds to 
Aigusy = 1212.9 GeV. It is obtained by varying the third generation scalar mass 
mo[I7J at the unification scale: The solid line is the formula, Eq. (|22p that is de- 
signed to numerically reproduce the effective Lagrangian result. The dashed lines are 
obtained from the formula (T23j) with inputs (m^ , ) = (m^ B + 5m, B + 5m) 

with various 5m values corresponding to M susy (= tl ~t, — — ) = 0.6,0.8, 1.0, 1.4 TeV. 
m^ = 125.5 ± 1 GeV is shown by blue band. 

natural SUSY, X t is approximately A t for A t ~ TeV. We should note that 
variations of the masses of the 1st and 2nd generations and gauginos from 
the NS3 inputs have little effect on since they are heavy in Natural SUSY 
scenario. 



1 The SOFTSUSY[82], SPheno[MIEU and SuSpect[85] codes use the same algo- 
rithm as Isajet[80 81] and employ similar threshold transitions matching the MSSM 
to the SM. The four codes produce mass spectrum in the mSUGRA model that 
are in close agreement. The Isajet [B0|81| code provides the NUHM2 model of our 
interest. 
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The rrih effective Lagrangian result with the NS3 input parameters can be 
numerically represented by the formula 



3mf 



m h = m hB (x t ) = M z c 2Pb + 2 



c + (ci + c 2 x t )X t 



X t = 2x 2 Al 
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M. 



(22) 



susy,B 



where the subscript B means the NS3 Benchmark point: C2p B = cos2/3# is 
calculated from tan/3 B = 19.4. M susy B is the SUSY breaking scale correspond- 
ing to (m t - 1)B ,m t - 2)jB ) 



(812.5, 1623.2) GeV; M susyiB 
1212.9 GeV. The coefficients 



(812.5 + 1623.2)/2 



(c 0)Cl;C2 ) = (2.661,0.2874,0.01717) 

have been determined by a least-squares fit with some weighting of the maxi- 
mal m>h region. 

We use m^B of Eq. (1221) as our benchmark at a given value of Xt- rrih values 
with different 2 and M susy = mtl ^ mt2 inputs are considered to be given with 
sufficient accuracy by shifting from rrih.B with a common value of x% through 
2- loop analytic formula ffl4l) . 



m 2 h (m h ,m i2 ,x t , tan/3) = 

m l,B^t) + [^h,EFT2 {m~tv m hi x ti tan P) ~ m h ! EFT 2 ( m i 1 ,B> m t 2 ,B, ^t, tan/3 B )] (23) 

In order to estimate the intrinsic uncertainty, we also consider the other nat- 
ural SUSY benchmark points, RNS1 and RNS2jT7], where rrih is estimated by 
using Isajet 7.83. 



mix 


m h 




A? s 


tan/3 


rrih (Isajet) 


RNS1 1416 


3425 


2420 


3764 


10 


123.7 


RNS2 1843 


4921 


3382 


5054 


8.55 


125.0 



m fe (Eq.((23 
124.1 
123.4 



(24) 



Here the masses and the Af s are given in units of GeV. The predictions from 
Eq. f[2"3"j) are given in the final column. Our formula (|2"B"|) is made by using 
a special input of NS3 benchmark point with M susy ~ 1.2 TeV, but it can 
be applied to wide range of cases with fairly good accuracy. The theoretical 
error of Eq. (|23|) is conservatively considered to be 2 GeV in whole range of 
parameters in natural SUSY scenario. 

In order to see the M susy dependence of m h , we shift the rn^ i2 from the NS3 
benchmark values commonly with 5m. The results are shown by dashed lines 



10 



in Fig. |2j which suggests the necessity of the maximal mixing condition when 
X t ~ v6M, 11 «, y |86J87| . The peak value of rrih gradually increases with ~lnM susy . 
The Higgs mass constraint rrih > 124.5 GeV requires a SUSY breaking scale 
M susy £ 0.6 TeV. 

The M susy dependence of m h in Natural SUSY points following Ref.[TT] are 
shown in Fig. [3j The points indicates a In M susy dependence, and in order to 
explain rrih > 124.5 GeV, it is indeed plausible that M susy > 1 TeV. 
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Fig. 3. M susy dependence of rrih i n Natural SUSY points following Ref.[T7]. The 
points are obtained from a scan over GUT scale parameters: the common scalar mass 
of the first two generations mo(l, 2) : 5 — 50 TeV, the third generation squark mass 
mo (3) : — 5 TeV, the common gaugino mass mu 2 : — 5 TeV, —4 < A t /mo(3) < 4, 
m A : 0.15 - 2 TeV, tan/3 : 1 - 60. See, \vL\'\j . 

The maximal mixing condition |Y 4 os | ~ 2M S , which corresponds to \X t \ ~ 
a/6M s in the DR or MS scheme, can be obtained [881175] by RGE running 
from the SUSY-GUT scale, as illustrated for Natural SUSY in Fig. HJ note 
that A t < is almost absent. The generated points are mainly in the region 
< At < 2; however, although improbable from the scan, the maximal mixing 
X t = v6M susy is possible in Natural SUSY. 

By taking rrih = 125.5 ± 2. GeV as a constraint to Eq. fl23|) . we can determine 
the allowed region in (TO£ 1 ,m£ 2 ) plane for a given value of 6~ t . Here we allow a 
somewhat large uncertainty of rrih,, 2 GeV, because of the theoretical uncer- 
tainty of our formula ( f23l) . The Higgs mass constraint severely constrains the 
top squark sector parameters, especially in that Am 4 -(= m t ~ 2 —m^) has a lower 
limit. From an Isajet scan over GUT scale parameters, we obtain the 9i de- 
pendence of Ami m Fig- El Almost all data points have large 0j, 1.3 < % < |, 
which means t\ ~ i R . Am t - decreases as % decreases from |. Actually 9i has a 
lower limit of 1.1 and we find that the on shell stop mass difference is bounded 
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Fig. 4. A t (M snsy ) / M susy dependence of in natural SUSY scan points. 
Xt = At — ucot(3 ~ At since /i is small, 150 GeV. The maximum of is not 
obtained at A t (M SVLSy ) / M SVLSy = V6 but at about 2, which is due to the difference of 
renormalization prescription of ISASUSY program, on-shell(OS) renormalization, 
and the EFT approach using the MS scheme. See, ref . [67] . 
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Concluding Remarks We have studied the implications for the scalar top 
sector of the recent Tevatron My/ measurements and the LHC and Tevatron in- 
dications of a 125 GeV Higgs boson. We utilized the H3m package to evaluate 
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mh through 3-loops in an effective Lagrangian approach with RGE evolution 
from the GUT scale. Natural SUSY was assumed, for which the third gener- 
ation scalar quarks are much lighter than the multi-TeV masses of squarks of 
the first two generations and the Higgsino mixing parameter \x is small, 150 
GeV. A maximal Higgs mass is attained that is close to the LHC experimental 
indications. The condition for maximal Higgs mass is an off-diagonal value of 
the stop mixing matrix X t = \/6M s in the DR renormalization scheme, which 
requires an on-shell soft-SUSY parameter at the weak scale of A t (M susy ) « 2 
TeV. The minimum value of the mass-splitting of two top squark states was 
found to be 400 GeV. As can be seen in Fig. [U the allowed region from the 
rrih constraint (blue region) satisfies the Mw constraint at 90% Confidence 
Level, independent of the value of % For 9i = ? a top-squark with sub- TeV 
mass is somewhat favored by the M w data; m$ < 500 GeV is possible for 
almost all values % when ti — £r. Precise experimental determination of 
at the LHC will tighten the restrictions on the top squark masses. The detec- 
tion of the scalar top states at the LHC would establish the SUSY theoretical 
underpinning of electroweak symmetry breaking. 
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